Quantum behaviour of a flux qubit coupled to a resonator 
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Q ■ We present a detailed theoretical analysis for a system of a superconducting flux qubit coupled to a trans- 

om ' mission line resonator. The master equation, accounting incoherent processes for a weakly populated resonator, 

is analytically solved. An electromagnetic wave transmission coefficient through the system, which provides a 
tool for probing dressed states of the qubit, is derived. We also consider a general case for the resonator with 
more than one photon population and compare the results with an experiment on the qubit-resonator system in 
the intermediate coupling regime, when the coupling energy is comparable with the qubit relaxation rate. 
-j_J ■ Keywords: superconducting qubit, transmission line, resonator. 
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C ' I. INTRODUCTION 

o : 

| Modern state of the art fabrication using nanotechnology brings together quantum optics and mesoscopic solid state physics. 
Different types of Josephson-junction quantum bits (qubits) - macroscopic quantum objects - are now intensively studied and 
their quantum behavior have been experimentally demonstrated (for review see e.g.01-0]). A series of quantum phenomena such 
as, for example, entanglement JHEjl], Rabi oscillations I7l- [l2tl . spin-echo and Ramsey fringes lfT3l[l4ll . Landau-Zener-Stiickelberg 

' interferometry fl5l4l9ll . have been recently demonstrated. Now, great interest is attracted to physics of artificial atoms (built on 
the basis of qubits) in a confined fields of electromagnetic resonators, which is known as circuit quantum electrodynamics 

• (CQED) i20ll . In pioneering CQED experiments, the artificial atom was electrostatically coupled to a high-quality transmission 

I line resonator. The large electrical dipole moment of the qubit and high energy density of the resonator allow this system to 
reach the strong coupling limit. This regime was studied theoretically II21H25I1 and experimentally for the charge qubit coupled 
capacitively to the resonator J26l - l28ll . Later inductive coupling for the flux qubit was proposed [29] and realized experimentally 

t-h ! ism. 

In this paper, we analyze the system of the superconducting flux qubit coupled to the transmission line resonator. Our aim is 
first to present the detailed theory of the qubit's states, dressed by the interaction with the quantum resonator, and their influence 
on the observable transmission. Second, we describe the regime of intermediate coupling studied recently experimentally by 
Oelsner et al. 13111 . Accordingly, the paper is organized as follows. In the next Section the model of the system is described. 
In Sec. Ill, we calculate energy levels of the system. Allowed transitions between the levels can be measured by spectroscopy 
applying external driving fields. Different representations of the system Hamiltonian are discussed in Sec. IV, particularly, 
the rotating-wave approximation (RWA), convenient for finding the stationary solutions. The analytical solution for the master 
equation is presented in Sec. V for the case of a weak and numeric calculations for strong driving regimes. In Appendix we 
present details of the theory for the transmission through the resonator. 
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II. DESCRIPTION OF THE SYSTEM 



We consider the flux qubit coupled inductively to a coplanar waveguide resonator, see Fig.Q] The flux qubit is a supercon- 
ducting loop with three Josephson junctions [32]. Two qubit states are naturally described in a flux basis. The two flux states 
(| t)> | 4-)) are differed by directions of the circulating current (clockwise and counterclockwise) in the loop. The qubit current 
interacts with the field of the resonator. The coplanar waveguide resonator is defined by two gaps in the transmission line, which 
form capacitances Co at x = ±1/2. The qubit is situated close to the center of the resonator (x = 0), where the current of the 
resonator fundamental mode is maximal. Note that the qubit dimensions are significantly smaller than the resonator wavelength, 
therefore we consider it as a point-like object. 
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FIG. 1: (a) Schematics of the qubit (denoted by a blue box) coupled to the transmission line resonator via inductance M. (b) Equivalent circuit 
for the description of the infinitesimal piece of length Ax of the transmission line, (c) Flux qubit with 3 Josephson junctions. 



The total Hamiltonian of the driven system 

H = -ffqb-r + 5(iw (1) 

is a sum of the driving field Hamiltonian H^ w and the qubit-resonator Hamiltonian 

-ffqb-r = -Hqb + H T + Hi ntl (2) 

which consists of the bare qubit i? q b, resonator H T and the interaction term i?j n t. The flux qubit Hamiltonian in the flux basis 
has the form JH 



#qb = --^ a x - 2°"*> ( 3 ) 

where A is the tunnelling amplitude, the energy bias e = 2/ p ($ — $o/2) is defined by the magnetic flux $, I p is the persistent 
current, a x , z are the Pauli matrices {a z \1) = — the current operator is 7 q b = —I p a z . 

The qubit is coupled to the transmission line resonator. A detailed resonator description is presented in Appendix (see also 
Refs. B3Ul33ll34ll ). The single-mode resonator is described by the following Hamiltonian 



H T = huj, (a^a + , (4) 

where a and are the annihilation and creation operators, which act at the number (Fock) states according to a \n) — y/n \n — 1) 
and \n — 1) = \pri \n). 

The term, describing the interaction between the resonator and the flux qubit, is 



H in t = MJ(0)J qb = -%(a f + a)a z , 



(5) 
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kg = MZro/p, (6) 

where M is the mutual loop-resonator inductance, 1(0) = I T o(a + (r) is the transmission line current operator (Eq. ( IA24b ). at 
the qubit's position, x — 0. 

The transmission line resonator is driven by the external probing voltage field at the frequency Wd close to the resonator 
characteristic frequency to x , as described in the Appendix. The qubit in turn is driven by the resonator photon field with the 
amplitude £ and the frequency o>d- The Hamiltonian of this field, described by photon exchange between the resonator and the 
driving field, can be written as 

= £ {iae 1 ^ 1 - iate-*"**) (7) 

(the derivation is presented in Sec. IV) 



III. ENERGY LEVELS AND THE SPECTROSCOPY OF DRESSED STATES 



The qubit-resonator Hamiltonian in the qubit eigenbasis (see, e.g., Ref. B36T1 ) can be written as 

#qb-r = H Q + H 'inti ( 8 ) 

where 

H' Q = ^a z +t^ T Ua+^\, (9) 



H' int = -%( a t + a ) (t^<J z ~ ^-a x ) , (10) 

with the bare qubit energy splitting 

hw qh = VA'+e 2 - (11) 
The bare system eigenstates are \e/g, n) = \e/g) (g) \n) and eigenvalues 



S e/g , n =±^ + ^ r ln+- 1 . (12) 

The states |e, n) and \g : n + 1) are degenerated at a; r = w q b and the degeneracy is lifted by the qubit-resonator interaction. The 
transition matrix element due to the interaction is 

(g,n+l \H[ Dt \e,n) = (e,n\H' int \g,n+ 1) = hg e Vn + l, (13) 

where the qubit-resonator interaction energy is 

A 

hg £ = hg- — . (14) 



Note that the coupling strength is scaled as A/fiw q b 112 ill . The eigenvectors |+) and |— ) of the total Hamiltonian H' qh l are 
obtained from the non-interacting qubit-resonator basis by the following transformation 



\ — ,n) \ _ f sinr] cosri ^ f \g,n + 1) 
|+, n) I \ — cos 77 sinrj 



(15) 



where 



2g £ Vn+l 

tan2?7= , (16) 



E± <n = fku r (n + 1) ± — 



(17) 



The energy of the ground state, \g, 0), is given by 



n n = V4g? (ra+l) + J 2 , 

6 = W q b — W r < 0. 

F -F M 
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(18) 
(19) 

(20) 



Here f2 n defines the energy difference £+. n — = Ml n . In particular, the energy anticrossing takes place at 8 = 0, that is 

at Suable*) = fiw r , and it is given by 



0.1 



A 



Q n (e*) = 2g e , v^Tl - 2g— VtTTT. 

nui r 



For example, in the inset of Fig.fSfa), the energy anticrossing is shown for n = 0. 
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FIG. 2: Energy diagrams calculated for the following qubit parameters A/h = 1.8 GHz, g/27r = 3 MHz, lu t /2it — 2.5 GHz. (a) Energy 
levels versus energy bias e. Avoided level crossing is shown as a close-up in the inset, (b) Contour lines of the energy difference versus bias e 
and the driving frequency oj^. The green (lower) line is for ftua = i5-,o — E sr and the black (upper) line is for ftuj^ — E+.o — E gI . (c) The 
same plot as in (b) but in a narrow vicinity to the resonator fundamental frequency w r . 



If the resonator is driven by a weak external field, so that it is weakly populated, one can limit the consideration by a few Fock 
states, neglecting unpopulated levels; the energy levels are plotted with Eqs. dTvT > and < f20b in Fig. [2ja). With the weak driving 
field the spectroscopy of the "dressed" energy levels can be done: the transmission is resonantly increased when the driving 
photon energy Hcj^ matches the system energy difference of Eqs. ( fT7b and d20l i, shown by two arrowed lines in Fig.[2fa) for two 
possible transitions. One can plot then the respective energy contour lines to describe experimental results (see Figs. |2jb, c)), 
which relates to the experimental data presented in Figs. 2 and 3 of Ref. HU]. With increasing driving amplitude, higher order 
processes such as multi -photon transitions may become possible l30l l37ll . 
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IV. HAMILTONIAN OF THE SYSTEM 

Jaynes-Cummings Hamiltonian. Let us rewrite the interaction Hamiltonian, Eq. ( [Tol l, by introducing the qubit lowering and 
raising operators 

a ± = i(cr x ±ia y ), (22) 

so that a + \g) = |e), a + |e) = 0, etc.; then we have 

H' int = hg e {a)<7~ + a<7 + ) + hg e {ao~ + a) ' a + ) - hg £ (a t + a)a z . (23) 

In vicinity of degeneracy of the states |e, n) and \g, n + 1), the second and the third term in Eq. ( 1231 can be neglected as they 
correspond to the processes, which require large extra energy. The first term together with H' from Eq. (O give the Jaynes- 
Cummings Hamiltonian 

ffjc = + fuvr (rfa + + Heifer' + aa+). (24) 

Interaction representation. We consider H' int in the interaction representation. For this, we note the following relations (see 
e.g. M) 

e MW ae- wW = ae~ iut , (25) 
e l t t ^ cr - e -^ t ' T2 = a~e~ iut . (26) 

Then we obtain 

H( n% = ei^H'^e-^ = hg e (aa + e^'^ + /i.e.) 
+ hge ( aa-e- i( ^ b+u ^ + h.c.) - hg-^— (ae' lu>A + h.c.) . (27) 

\ / huj qh 

In the RWA, when w q b — Lu r <C w q b, the first term is slowly rotating, while the second and third terms are fast rotating ones. 
This justifies neglecting these terms. 

Driving Hamiltoinan. We will consider scattering of the right propagating wave y{ (e~ lk ^- x+l ^ 2 ' >+ ' lu ' dt + c.c.) on the resonator 
(see Eq. ( 1A31| >), where V{ is chosen to be a real amplitude. This wave drives the resonator, which in turn generates the scattered 
waves. Using semiclassical approach, we first consider the driving dynamics under the classical field and then calculate the 
field generated by the resonator. The calculated first order scattering gives an exact solution because the two scattered waves, 
propagating in different directions cancel out the second order driving effect (see Eq. dA641 >). The dipole-like interaction 
Hamiltonian can be presented as a product of the incident wave voltage field and charges generated by the resonator field on the 
coupling capacitances CoV T o(ia — ia)) sin (±k r l/2) (see Eq. jA27\ ) 

= { (e Mdt + e-^ dt ) (ia - tot) ; (28 ) 

where 

£ = C V{V M sin (M/2)(-l + e-' iM ) = -2C Q V{V r0 . (29) 

And omitting fast rotating terms in RWA, we arrive to Eq. ((7). In these equations we assume B\ = ioCqZi <C 1, which is valid 
for high quality resonators. 

Rotating-wave approximation. We consider the Hamiltonian of the driven system in the RWA 

#rwa = U (H' qb _ T + H^) C/t + ihUtf . (30) 

For this we choose the transformation 

U = exp [iw d t (a f a + a z /2)] (31) 
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and obtain 

Hrwa — h — 2~° z ^ LU ^ a ^ a + ^Se (acr + + o'er - ) + £(ia — iaJ), (32) 

Suj qh = W qb - UJd, (33) 

<5u r = w r — Wd- 

Control microwave field. For the sake of generality, we consider also the case when the qubit is driven by the separate 
microwave field, coupled, for example via an additional microwave line. Then, we have 

= - J P $ ac cos OJ d t ■ a z , (34) 

where $ ac is the amplitude of the driving flux. In the qubit eigenstate representation, this is simplified to the form of Eq. (0 
according to 

iu) A t i p ~iu d t f - A \ 

«ff = (j- *. - - {.(e-V- + e-V + ), ,35, 

Here we have left only slowly rotating terms (see discussion above). Note that the amplitude £ e is dependent on the bias e (see 
Eq. (HU). Then in the RWA after the transformation d3~Tl i we obtain the expression which differs from Eq. d32l by substituting 
the last term with £, £ (<j + + <r~). 

Dispersive regime. In the dispersive regime (when 5 ^> hg e ) the diagonalization of the Hamiltonian ( 1241 in the second order 
in g/S 113 811 gives 

H = ~ (ftu qh + <?z + ffi^r + ~^ a ^j a^a. (37) 

This expression explicitly shows the qubit transition energy shift by the coupling and also the resonator energy shift by the qubit, 
which sign depends on the qubit state. 

V. SOLUTION OF THE MASTER EQUATION FOR THE DENSITY MATRIX OF THE SYSTEM 

To describe the qubit-resonator dissipative and incoherent dynamics we assume that all processes in our system are Markovian 
and solve the master equation for the density matrix p 

p = - t -[H,p]+C[p}. (38) 



It includes the dynamic part and dissipative Lindblad term [ 3S 

, 3 



k=l 

where 



Ci = Sn*, 7i = ^, (40) 

r - P+n n, - 1 - 1 1 
V I lcf> -1 2 *-Ll 

C3 = \J~Ha. 
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The Lindblad operator C presents dissipation in the resonator (photon decay) with the rate x = x ex t + XinU where x ex t and x- mt 
are external (leaking out through of the resonator) and internal (resistive) loss rates, and the qubit decoherence consisting of the 
relaxation rate 71 and the dephasing rate 7^. Here we consider nondispersive regime (near the qubit-resonator resonance). The 
Hamiltonian of the system H in the rotating wave approximation has the form of Eq. (f32t . The solution of the master equation 
determines the observable quantities, in particular, the expectation value of the photon field in the resonator 

(a) = Tr(ap). (41) 

The Hilbert space of the composite system is the tensor product of the qubit space and the photon space with basis vectors 
\e/g,n) — \e/g) <8> \n). Basis vectors \g) and |e) 
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are the eigenvectors of the operator a z . Vectors of Fock (photon) states \n) (the eigenvectors of the photon number operator 
a) a \n) = n \n)) are the vectors in the infinite-dimensional space N — 00 
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In the basis \e/g, n) the matrix equation d38l ) is the infinite set of equations for the infinite-dimensional matrix py. 

Below, we consider the simplest case of N = 2 (weak driving limit), where the analytical solution is possible, and in the case 
of N ^ 1 we study the problem numerically. 



A. Weak driving limit 

To find the analytical solution we restrict the photon space to N = 2, assuming that the mean photon number in the resonator 
(created by the driving field with the amplitude £) is much less than unity. The basis \e/g,n) in this case consists of 4 base 
vectors hi 

&! = |.g0) , & 2 = |e0) , 63 = \g\) , b 4 - |el> , (44) 
and the density matrix = (bi \p\ bj) takes the form 



/ PeCUO PeO,gO PeO,el PeO,gl 
_ PgO,eO PgO.gO PgO,el PgO.el 
Pel,eO Pel.eO Pel, el Pel,gl 
V Pgl,eO Pgl.gO Pgl,el Pgl.gl 

In the steady state from Eq. d38l > we have 16 linear equations for the matrix elements py. In the weak driving limit, when the 
drive does not change population of the ground state p g o, g o = 1, leaving up to the first order terms only in the amplitude £, we 
obtain the density matrix p. The nonzero elements of the matrix p^ are 




PgO,gO - 1, 

_ * -im)se 

Pe0,g0 ~ P 90 ,e0 " g 2 _ ^ _ . f )( ^ b _ , 



where 7 = + 7^ 
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Using d45b in fiTt we obtain for the mean value of the voltage field in the resonator in the weak driving (WD) limit for positive 
frequencies 



£(<!>Wqb + Vy) 



g 2 e - (Su t + if)(Su} qh + i-y) ' 



(47) 



The transmission coefficient of the output driving signal t is defined by the photon field in the resonator, Eq. ( IA63I ), and according 
with (|47| | we obtain 



*WD- 



2 gf-Oftjr + if)(5u qh +i~f)' 



(48) 



When g £ = 0, this equation gives the transmission coefficient through the resonator, which for the linear resonator can be derived 
classically (Eq. |A40V The plot of the transmission amplitude |t| WD , given by Eq. d48| i, is shown in Fig.|3]for u; q b = ui Y and 
different values of the decay rates k and 7 (given in units of the coupling constant g E ). For weak decay rates yt and 7, the 
transmission spectrum displays the Rabi-splitting peaks (red solid curve), which are smeared with increasing of the decay. 
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FIG. 3: Normalized transmission amplitude \t\ as a function of the driving frequency detuning — u; r at e = e* (when tj q b(e*) = to T ) for 
different values of the decay rates x and 7 (given in the figure in units of g e ), calculated with Eq. J48t . 
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FIG. 4: Normalized transmission amplitude \t \ as a function of the bias e and the driving frequency detuning Ud — Ut, calculated with Eq. d48t 



In Fig. @] the density plot of the transmission amplitude as a function of the bias e and the detuning uja — cj t is shown. The 
parameters here and below are taken for the comparison with the relevant experimental work 13111 A/h = 1.8 GHz, g/27r = 3 



9 



MHz, u> x /2tt — 2.5 GHz (the same as in Fig. [2j and also the loss rate of the resonator x/2tt = 1.25 ■ 10~ 4 GHz and the loss 
rate of the qubit 7 = g. Note that we consider the intermediate coupling regime, when g = 7 ^> x. The transmission amplitude 
is resonantly increased along the lines shown in Fig. |2c) as expected. In the narrow vicinity of the resonator characteristic 
frequency, u>d € (w r — g, u r + g), the avoided crossing at e = e* is demonstrated, as it was reported in Ref. Bill . 

For more detailed comparison and finding the parameters with better accuracy (e.g. decay rate 7), we need to compare 
experimental and theoretical sets of crossection of surfaces \t\ versus e and u>d- This is shown in Fig.[5]for uid = to T - 



theory 

experiment 




I 1 1 1 1 1 1 1 1 

-4-3-2-101234 

e 



FIG. 5: Normalized transmission amplitude \t\ as a function of the bias e for tuj = oj t calculated with Eq. d48l >, red line, and obtained 
experimentally, blue line. 



B. Numerical solution of the master equation. Beyond the weak driving regime. 

In the case, when driving is not weak, i.e. the mean photon number (a^a) > 1, we have solved the equation for the density 
matrix p numerically. The results are presented in Fig. [6] The transmission amplitude \t\ in all cases is normalized on the 
maximal value at w q b = ui T . In Fig. [(J a) the transmission amplitude is shown for the case of small damping x/g £ = 0.1 and 
7/ge = 0.1. At a weak driving amplitude £ the red curve in Fig. [6ta) coincides with \t\ WD (cud) (Fig. 01. With increasing £, each 
split Rabi peak is additionally split (blue curve) (see also in Ref. L40j). With further increasing of the amplitude £, the additional 
splitting is smeared (green curve). Thus in the nonlinear regime we observe the qualitatively new features as compared to the 
weak driving limit. When the decay is rather large, such that in the weak-driving case, we do not have the Rabi splitting (green 
curve in Fig. |3), in the nonlinear response, we do not observe the qualitatively new features, as shown in Fig. [6jb) (x/g E = 1 
and 7/g e = 2). 

We also calculate the average number of photons in the resonator, n = (a^a). For the parameters in Fig. [6] it depends on the 
frequency; the maximal values are the following n max = 0.005 for £/g e ~ 0.01, n max = 0.3 for £/g e = 0.15, n max = 1-8 for 

e/g £ = 0.3. 

VI. CONCLUSION 

We presented the detailed theory for the system of the flux qubit coupled inductively to the transmission line resonator. The 
transmission coefficient is calculated with the system's density matrix by solving the master equation within RWA. 

The avoided crossing of the dressed energy levels is shown in the resonant case, where ujd ~ ^>ab ~ w r . This is demonstrated 
in the intermediate coupling regime, which describe the experimental results of Oelsner et al. Olll . We have shown that the 
dissipation smears the Rabi splitting. Moreover, we have demonstrated the double splitting in the strong driving regime. 
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Appendix A: Transmission line resonator 



In this Appendix we consider the resonator formed by the transmission line interrupted by two capacitances Co- The qubit we 
assume to be coupled inductively to the resonator at its center, see Fig. [Ha). We start by presenting the equations which describe 
the transmission line. 



1. The transmission line 



The transmission line is usually modelled as an infinite series of the elementary circuits (e.g., 14111 ). as shown in Fig. \V[b). 
Here elementary inductance, capacitance and conductance are AL = LAx, AC = CAx, AG — GAx, where L, C and G 
are inductance, capacitance and conductance (of parallel resistance) per unit length. For the circuit in Fig. [TJb), we can write 
(neglecting the Ohmic losses) the equations for the transmission line, by applying the Kirchhoff 's laws for the voltage V(x,t) 
and the current I(x, t); in the limit Ax — > they take the form 

dV(x,t) = L dl(x,t) 
Ox Ot 

These equations can be rewritten for either I(x, t) or V(x, t) as following 

2 A 1 8 2 A _ xdA 

dx^~^^-^~8t> A-{J,V}, (A3) 

(A4) 
(A5) 

Here v is the phase velocity and k defines the loss in the transmission line. 
Assuming I(x, t) — I(x)e iuJt and V(x, t) = V(x)e lujt , we obtain 




dV(x) 

dx 
dl(x) 
dx 



-iu)LI(x), (A6) 
-(G + iuC)V(x). (A7) 
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Then equation for A(x) = {I(x), V(x)} can be written as following 



££W-j>A{x) = 0, (A8) 



7 = y/iu)L(G + iuC) = a + ik. (A9) 
Solving equation for V(x) and using Eq. (IA6t . we obtain 

V(x) = V^e-^ + V Q l e^ x , (A10) 

I(x) = ^e-^-^-e^, (All) 



where 



Z = — = Z 1 +iZ 2 . (A12) 
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* = ^ = (M3) 



fc 2 ' 



When losses in the line are small (G <C wC), we obtain 



knwVW=-,an—\— = £-, (A14) 



2 V C 2u' 



Zi = \l^, Z 2 = —. (A15) 



Here the constants V r and Vq are the amplitudes of the right- and left-moving waves and Zq is the transmission line characteristic 
(wave) impedance. 

2. Open transmission-line resonator 

We consider the open transmission line of the length I. The quality factor of the resonator 114 111 can be written as 

„ k UJ r C L0 r 

Q = 7r = -$r = —. (A16) 

2a G x 

Now let us define normal modes of the resonator without dissipation (x = 0). Then assuming zero current through the boundaries 
at x = ±1/2 for this modes, we obtain 

lU)(x) = (e-^ x - (-l)V* a ) , (A17) 
Zo 

V&(x) = V r (e~ ik ^ + {-lYe***) , (A18) 
where kjl = jir, j = 1, 2, 3, .... In particular, for the fundamental mode j = 1 of the resonator we obtain 

2V r 

I (1) (x) = — -coshx, (A19) 
Zq 

V {1) (x) = -2iV r sin kix. (A20) 
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For the fundamental mode j = 1 of the A/2 resonator (I = A/2), we have fc r = fci = 7r/Z, cj r = wi = fciv = ^-^=, where 
L r = LI and C r — CI are the total inductance and capacitance of the resonator. 

Quantization of the resonator eigenmodes results in the following expressions for the current and voltage operators and the 
Hamiltonian 



We consider the frequency close to the fundamental mode frequency u; r , and, therefore, we ignore other modes. For the 
fundamental mode, with fci = ir/l and omitting the index j — 1, we obtain 



I = I^a + at) cos™ , J r0 = J^, (A24) 



V = l Ko(a-a t )sin^, V l0 = ^^, (A25) 

where I T o and V x o are the zero-point root mean square (rms) current and voltage, and the Hamiltonian is given by Eq. |@). 

We also consider the realistic case: the resonator with two point-like coupling capacitances Cq at the ends with 6i — 
ujCqZi <g; 1. For the fundamental mode, the current and voltage operators are modified to 

/ = Iro(a + a T ) cosfc r x, (A26) 
V = V T o (ia — ia t) sin k T x. (A27) 

And we find 

fc r w fci - — i (A28) 

from the boundary conditions 

I l0 (a + a f ) cos(±fc r ;/2) = K. [ia(-«wC ) - m^wCy sin(±fc r Z/2). (A29) 
This results in the shifted resonant frequency 

( 26>! \ 

Wr = wi 1 , (A30) 



7T 

which is slightly lower than the fundamental frequency due to external coupling to the outside lines via the capacitance Cq. 



3. Transmission through the coplanar waveguide resonator 

Now we will consider a classical problem of transmission of waves through the resonator. It will help us to find correspondence 
between the classical and quantum-mechanical solutions and to define the photon decay rates. The incident wave propagates 
from left to right and interacts with the transmission-line resonator at a; = — 1/2 through the capacitance Cq. The output wave 
is detected after another capacitance Cq at x = 1/2. We will obtain the system of equations for VJ and Vj, which define the 
classical current and voltage in j-th region, j — 1, 2, 3, respectively for x < — 1/2, x e {—1/2, 1/2), and x > 1/2; 

Vj(x) = VTe-rt*-*** + V]ei {x - X >\ (A31) 
V r V 1 

Ij(x) = -Z-e-* 1 ^-**) - -J-e 1 ^-^, (A32) 
Zq Zq 



13 



where x\ — —1/2, x 2 = and xj, = 1/2, We assume the matched termination (with impedance equal to Z$), then there is no 
left-propagating wave in the third region, V 3 ' = 0. The boundary conditions for currents and voltages at the points x = ±1/2 are 
the following 

h(-l/2) = I 2 (-l/2), (A33) 
I 2 (l/2) = / 3 (Z/2), (A34) 



V x {-l/2) = V 2 (-l/2) + I 2 {-l/2)/iuC , (A35) 

V 2 (l/2) = V 3 (l/2) + I 3 (l/2)/iu>C . (A36) 

From these equations, substituting Vi(— 1/2) = V{ + Vj, we find a useful relation between the field in the resonator and the 
external field V3 

V 3 (l/2) = V 2 (l/2)-^-, (A37) 
1 + 1&1 

where 6\ = ujCqZ^. 

We define the transmission coefficient t as a ratio between the transmitted wave and the incident one as 

yr 

t = ^ (A38) 

and find directly from Eqs. dA3UA35b 

40?e"T ; 

t ~ 46*2 - Ai6 x - 1 + e-W ■ (A39) 

For the interesting case of high-Q resonators (a <C k and 9 ss loCqZ\ <C 1), we can express the transmission coefficient near 
the fundamental mode (uji = vn/l) in the compact form 

^rfw' (A40) 



where 5lu = w T — lj^ is detuning from the resonant frequency 



Wr = Wi ( 1 - — ] (A41) 



7T 

due to coupling capacitance Co- Note that this formula coincides with Eq. ( IA30b . however, obtained from the classical solution. 
The peak width 

Acj = x, (A42) 
is determined by the total photon decay rate x = xr cx t + Xint which is the sum of the photon decay rate due to the external loss 

_ 46^w r 

Xoxt — 

7T 

determined by the coupling to the external transmission lines via Co and the internal photon decay rate 

2alcu T 



(A43) 



due to dissipations within the resonator. The quality factor is 

: s 



(A44) 



Q = ^ = WT^i- (A45) 



This rate is consistent with its definition given in Ref. 13 111 . 

Below we estimate the photon decay rate x for the coplanar waveguide resonator with parameters taken from [42] / = 23 
mm, w r /27r = 2.5 GHz, Co = lfF, Z\ = 50 Ohm, which give B\ = 7.8 x 10~ 4 . The capacitance per unit length C is calculated 
from the expression for B\ at resonance Q\ = ttCq/IC. We thus obtain C = 1.74 x 10~ 10 F/m. Finally, for the photon decay 
rate x we obtain >c/2tt = 1.95 kHz. This value is about two times smaller than the ones obtained in [42]. We assume that this 
discrepancy is due to dielectric losses G. It allows us to estimate a from AB\ « 2al, then a ~ 5.3 x 10~ 5 m _1 . Therefore, for 
G we obtain G = 2a/Z 1 sa 2.12 Ohm" 1 m" 1 . 
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a. Transmission in the dispersive regime 

Here we consider an effect of the qubit on the transmission coefficient, substituting the qubit by an additional classical 
inductance coupled to the resonator. This classical analogy may be helpful to understand the quantum-mechanical effect. In the 
dispersive regime, coupling to the qubit can be described as an additional classical inductance L q b at the position x — 0. Such 
a problem is described by adding two more equations for x = to the system of equations ( lA33HA36l l. which follows from 
Eq. dAU by adding to the r.h.s. the following term 



xi \nr dI ^ x , , T dl{x,t) 
6(x)M—^- = -d{x)L qh - 



where 



In the ground state we have [43, 4- 



->qb 



Of 



M 



9$ ■ 



4M 2 /2A 2 

Lqb = 



( A 2 +£ 2)3/2- 



We modify the definition of V 2 given in Eq. ( [A3 It 



(A46) 



(A47) 



(A48) 



V 2l (x) = V 2 r l e-^ + Vi l e^ 



< x < 0, 
/ 



V 2r {x) = V 2 r r e-~< x + V\ r &* , < x < -, 



(A49) 
(A50) 



then the boundary conditions at x = are 

Jaj(0)=/ar(0), 

V 2 i(0)^V 2r {0) + iu J L qh I 2r (0). 
The solution of the system of equations for the transmission coefficient can be written as 

n -1 



t' = 



- l-i20i)' 



1 - i- 



Suj - 



}LqbK cx t \ 



(A51) 
(A52) 



(A53) 



where t is the transmission without the qubit (L q = 0) from Eq. ( |A39t . Here we used the following simplifications —7/ 
—ikl — al = —in — al and al -C 1, 6\ <€. 1. Finally we rewrite t 1 in the compact form 



t' 



1 - i 



2Su 



(A54) 



where detuning Soj' = uj[ — uj^ from the redefined resonance frequency 



(A55) 



which is shifted due to the extra inductance L qb in the resonator. The phase shift of the transmission coefficient t' at Soj = is 
found as 



tamp = 



Im[t'] _ w r L qb >c cxt _ 1 / C T \ L qb x cxt 



Re[t'\ 26\Z x 2tt \CqJ L r x 



In the ground state we obtain 



tan ip 



A 



l + (e/A)' 



-3/2 



(A56) 



(A57) 
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b. Resonant transmission 



The measured resonator field is expressed via the field operator expectation values (I) or (V) (see Eqs. (IA24I IA25t >). Par- 
ticularly, the expectation value of the voltage at x = ±1/2 is V t q (ia — ia' ) sin (±M/2) ~ ±V l0 ( ia — ia)\, and the positive 
frequency component of the charge on the capacitances Co are 

(<?+) =±C V T0 (-i^). (A59) 

The current leaking out from the resonator, expressed via the reflection and transmission coefficients r and t, is a time-derivative 
of the charge at x = —1/2 and x = 1/2, that is 

^L(l-r)=iw<g+). (A60) 

^t = ioj(q + ). (A61) 
"1 



Then the transmission coefficient can be presented as 



t = ^, (A62) 
which after some algebra using Eq. ( f29T > and Eq. ( IA44I ) can be rewritten in a simple physical form 

' = -#)<"'>• <A«) 

It is also straightforward to demonstrate that the scattered waves: back scattered V{re lk ^ x+l ^ 2 ^ at x < —1/2 and forward 
y^te~ lk{ - x ^ l l 2 ^ — V{e~ ik ( x ~ l ' 2 > (difference between the transmitted and the undisturbed one as, if there is no resonator) at 
x > 1/2 are equal in amplitude (1 — r = t) and, therefore, effectively result in zero interaction energy 

C V{V l0 sm{k r l/2)[-r + (1 - t)e~ ikrl ] =0, (A64) 

(compare with Eq. d28l )) that is in the quasi-classical approach of scattering, the first order scattering gives an exact solution. 
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